Rapid Decline of SARS-CoV-2 Viral Load in Single vs. Double-Dose (Short-Interval <6 Weeks) ChAdOx nCoV-19 Vaccinated Health-Care Workers.

The present work was carried out during the emergence of Delta Variant of Concern (VoC) and aimed to study the change in SARS CoV-2 viral load in Covishield vaccinated asymptomatic/mildly symptomatic health-care workers (HCWs) to find out the optimum isolation period. The SARS CoV-2 viral load was carried out in sequential samples of 55 eligible HCWs which included unvaccinated (UnV; n = 11), single-dose vaccinated (SDV, n = 20) and double-dose vaccinated [DDV, n = 24; short-interval (<6 weeks)] subjects. The mean load of envelope (E) gene on day 5 in SDV [0.42 × 105 copies/reaction] was significantly lower as compared to DDV [6.3 × 105 copies/reaction, P = 0.005] and UnV [6.6 × 105 copies/reaction, P = 0.001] groups. The rate of decline of SARS CoV-2 viral load in the initial 5 days of PCR positivity was significantly higher in SDV as compared to that in DDV (Mean log decline 0.39 vs. 0.19; P < 0.001). This was possibly due to interference of adenoviral immunity of first dose of adenovirus-vectored vaccine in double-dose vaccinated HCWs who had received vaccines within a shorter interval (<6 weeks).

Draft genome sequence of Kluyvera ascorbata HAK22 isolated from human feces.

The whole genome sequence of rare human pathogen Kluyvera ascorbata strain HAK22 is reported. The K. ascorbata HAK22 was isolated from healthy human from Gurugram, Haryana, India. The draft genome has a length of 4.7 Mbp, with 54.36% GC content and 4,411 proteins, 4,470 genes, and 18 antimicrobial resistance genes.

Computational study of the motor neuron protein KIF5A to identify nsSNPs, bioactive compounds, and its key regulators.

Introduction: Kinesin family member 5A (KIF5A) is a motor neuron protein expressed in neurons and involved in anterograde transportation of organelles, proteins, and RNA. Variations in the KIF5A gene that interfere with axonal transport have emerged as a distinguishing feature in several neurodegenerative disorders, including hereditary spastic paraplegia (HSP10), Charcot-Marie-Tooth disease type 2 (CMT2), and Amyotrophic Lateral Sclerosis (ALS). Methods: In this study, we implemented a computational structural and systems biology approach to uncover the role of KIF5A in ALS. Using the computational structural biology method, we explored the role of non-synonymous Single Nucleotide Polymorphism (nsSNPs) in KIF5A. Further, to identify the potential inhibitory molecule against the highly destabilizing structure variant, we docked 24 plant-derived phytochemicals involved in ALS. Results: We found KIF5AS291F variant showed the most structure destabilizing behavior and the phytocompound "epigallocatechin gallate" showed the highest binding affinity (-9.0 Kcal/mol) as compared to wild KIF5A (-8.4 Kcal/mol). Further, with the systems biology approach, we constructed the KIF5A protein-protein interaction (PPI) network to identify the associated Kinesin Families (KIFs) proteins, modules, and their function. We also constructed a transcriptional and post-transcriptional regulatory network of KIF5A. With the network topological parameters of PPIN (Degree, Bottleneck, Closeness, and MNC) using CytoHubba and computational knock-out experiment using Network Analyzer, we found KIF1A, 5B, and 5C were the significant proteins. The functional modules were highly enriched with microtubule motor activity, chemical synaptic transmission in neurons, GTP binding, and GABA receptor activity. In regulatory network analysis, we found KIF5A post-transcriptionally down-regulated by miR-107 which is further transcriptionally up-regulated by four TFs (HIF1A, PPARA, SREBF1, and TP53) and down-regulated by three TFs (ZEB1, ZEB2, and LIN28A). Discussion: We concluded our study by finding a crucial variant of KIF5A and its potential therapeutic target (epigallocatechin gallate) and KIF5A associated significant genes with important regulators which could decrypt the novel therapeutics in ALS and other neurodegenerative diseases.

Genome sequencing of SARS-CoV-2 omicron variants in Delhi reveals alterations in immunogenic regions in spike glycoprotein.

The SARS-CoV-2 omicron variants keep accumulating a large number of mutations in the spike (S) protein, which contributes to greater transmissibility and a rapid rise to dominance within populations. The identification of mutations and their affinity to the cellular angiotensin-converting enzyme-2 (ACE-2) receptor and immune evasion in the Delhi NCR region was under-acknowledged. The study identifies some mutations (Y505 reversion, G339H, and R346T/N) in genomes from Delhi, India, and their probable implications for altering the immune response and binding affinity for ACE-2. The spike mutations have influenced the neutralizing activity of antibodies against the omicron variant, which shows partial immune escape. However, researchers are currently exploring various mitigation strategies to tackle the potential decline in efficacy or effectiveness against existing and future variants of SARS-CoV-2. These strategies include modifying vaccines to target specific variants, such as the omicron variant, developing multivalent vaccine formulations, and exploring alternative delivery methods. To address this, it is also necessary to understand the impact of these mutations from a different perspective, especially in terms of alterations in antigenic determinants. In this study, we have done whole genome sequencing (WGS) of SARS-CoV-2 in COVID-19 samples from Delhi, NCR, and analyzed the spike's mutation with an emphasis on antigenic alterations. The impact of mutation in terms of epitope formation, loss/gain of efficiency, and interaction of epitopes with antibodies has been studied. Some of the mutations or variant genomes seem to be the progenitors of the upcoming variants in India. Our analyses suggested that weakening interactions with antibodies may lead to immune resistance in the circulating genomes.

Investigating the link between miR-34a-5p and TLR6 signaling in sepsis-induced ARDS.

Acute lung injury (ALI)/acute respiratory distress syndrome (ARDS) are lung complications diagnosed by impaired gaseous exchanges leading to mortality. From the diverse etiologies, sepsis is a prominent contributor to ALI/ARDS. In the present study, we retrieved sepsis-induced ARDS mRNA expression profile and identified 883 differentially expressed genes (DEGs). Next, we established an ARDS-specific weighted gene co-expression network (WGCN) and picked the blue module as our hub module based on highly correlated network properties. Later we subjected all hub module DEGs to form an ARDS-specific 3-node feed-forward loop (FFL) whose highest-order subnetwork motif revealed one TF (STAT6), one miRNA (miR-34a-5p), and one mRNA (TLR6). Thereafter, we screened a natural product library and identified three lead molecules that showed promising binding affinity against TLR6. We then performed molecular dynamics simulations to evaluate the stability and binding free energy of the TLR6-lead molecule complexes. Our results suggest these lead molecules may be potential therapeutic candidates for treating sepsis-induced ALI/ARDS. In-silico studies on clinical datasets for sepsis-induced ARDS indicate a possible positive interaction between miR-34a and TLR6 and an antagonizing effect on STAT6 to promote inflammation. Also, the translational study on septic mice lungs by IHC staining reveals a hike in the expression of TLR6. We report here that miR-34a actively augments the effect of sepsis on lung epithelial cell apoptosis. This study suggests that miR-34a promotes TLR6 to heighten inflammation in sepsis-induced ALI/ARDS. Supplementary Information: The online version contains supplementary material available at 10.1007/s13205-023-03700-1.

Betasatellite-encoded ßC1 protein regulates helper virus accumulation by interfering with the ATP hydrolysis activity of geminivirus-encoded replication initiator protein.

Geminivirus-betasatellite disease complexes are an epidemic threat to the majority of economically important crops across the world. Plant virus satellites including betasatellites are maintained by their associated helper virus. Geminivirus-betasatellites influence viral pathogenesis by substantially increasing or decreasing their helper virus accumulation. In the present study, we attempted to understand the mechanistic details of the geminivirus-betasatellite interaction. Here, we used tomato leaf curl Gujarat virus (ToLCGV) and tomato leaf curl Patna betasatellite (ToLCPaB) as a model system. This study reveals that ToLCGV can efficiently trans-replicate ToLCPaB in Nicotiana benthamiana plants, but ToLCPaB greatly reduced the accumulation of its helper virus DNA. For the first time, we have identified that the ToLCPaB-encoded ßC1 protein is able to interact with ToLCGV-encoded replication initiator protein (Rep). In addition, we demonstrate that the C-terminal region of ßC1 interacts with the C-terminus of Rep (RepC) protein. Our previous study had established that ßC1 proteins encoded by diverse betasatellites possess a novel ATP hydrolysis activity and the conserved lysine/arginine residues at positions 49 and 91 are necessary for this function. Here, we show that mutating lysine at positions 49 to alanine of ßC1 (ßC1K49A) protein did not affect its ability to interact with RepC protein. Biochemical studies performed with ATP hydrolysis activity-deficient K49A mutated ßC1 (ßC1K49A) and RepC proteins revealed that Rep-ßC1 interaction interferes with the ATP hydrolysis activity of Rep protein. Further, we demonstrate that ßC1 protein is able to interact with D227A and D289A mutated RepC proteins but not with D262A, K272A or D286A mutated RepC proteins, suggesting that the ßC1-interacting region of Rep protein encompasses its Walker-B and B' motifs. The results of docking studies supported that the ßC1-interacting region of Rep protein encompasses its motifs associated with ATP binding and ATP hydrolysis activities. Docking studies also provided evidence that the Rep-ßC1 interaction interferes with the ATP binding activity of Rep protein. Together, our findings suggest that ßC1 protein regulates helper virus accumulation by interfering with the ATP hydrolysis activity of helper virus Rep protein.

Monkeypox virus: phylogenomics, host-pathogen interactome and mutational cascade.

While the world is still recovering from the Covid-19 pandemic, monkeypox virus (MPXV) awaits to cause another global outbreak as a challenge to all of mankind. However, the Covid-19 pandemic has taught us a lesson to speed up the pace of viral genomic research for the implementation of preventive and treatment strategies. One of the important aspects of MPXV that needs immediate insight is its evolutionary lineage based on genomic studies. Utilizing high-quality isolates from the GISAID (Global Initiative on Sharing All Influenza Data) database, primarily sourced from Europe and North America, we employed a SNP-based whole-genome phylogeny method and identified four major clusters among 628 MPXV isolates. Our findings indicate a distinct evolutionary lineage for the first MPXV isolate, and a complex epidemiology and evolution of MPXV strains across various countries. Further analysis of the host-pathogen interaction network revealed key viral proteins, such as E3, SPI-2, K7 and CrmB, that play a significant role in regulating the network and inhibiting the host's cellular innate immune system. Our structural analysis of proteins E3 and CrmB revealed potential disruption of stability due to certain mutations. While this study identified a large number of mutations within the new outbreak clade, it also reflected that we need to move fast with the genomic analysis of newly detected strains from around the world to develop better prevention and treatment methods.

Combinatorial Network of Transcriptional and miRNA Regulation in Colorectal Cancer.

Colorectal cancer is one of the leading causes of cancer-associated mortality across the worldwide. One of the major challenges in colorectal cancer is the understanding of the regulatory mechanisms of biological molecules. In this study, we aimed to identify novel key molecules in colorectal cancer by using a computational systems biology approach. We constructed the colorectal protein-protein interaction network which followed hierarchical scale-free nature. We identified TP53, CTNBB1, AKT1, EGFR, HRAS, JUN, RHOA, and EGF as bottleneck-hubs. The HRAS showed the largest interacting strength with functional subnetworks, having strong correlation with protein phosphorylation, kinase activity, signal transduction, and apoptotic processes. Furthermore, we constructed the bottleneck-hubs' regulatory networks with their transcriptional (transcription factor) and post-transcriptional (miRNAs) regulators, which exhibited the important key regulators. We observed miR-429, miR-622, and miR-133b and transcription factors (EZH2, HDAC1, HDAC4, AR, NFKB1, and KLF4) regulates four bottleneck-hubs (TP53, JUN, AKT1 and EGFR) at the motif level. In future, biochemical investigation of the observed key regulators could provide further understanding about their role in the pathophysiology of colorectal cancer.

Amyotrophic Lateral Sclerosis Risk Genes and Suppressor.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease that leads to death by progressive paralysis and respiratory failure within 2-4 years of onset. About 90-95% of ALS cases are sporadic (sALS), and 5-10% are inherited through family (fALS). Though the mechanisms of the disease are still poorly understood, so far, approximately 40 genes have been reported as ALS causative genes. The mutations in some crucial genes, like SOD1, C9ORF72, FUS, and TDP-43, are majorly associated with ALS, resulting in ROS-associated oxidative stress, excitotoxicity, protein aggregation, altered RNA processing, axonal and vesicular trafficking dysregulation, and mitochondrial dysfunction. Recent studies show that dysfunctional cellular pathways get restored as a result of the repair of a single pathway in ALS. In this review article, our aim is to identify putative targets for therapeutic development and the importance of a single suppressor to reduce multiple symptoms by focusing on important mutations and the phenotypic suppressors of dysfunctional cellular pathways in crucial genes as reported by other studies.

Hereditary Vitamin D-Resistant Rickets (HVDRR) associated SNP variants of vitamin D receptor exhibit malfunctioning at multiple levels.

Vitamin D receptor (VDR) is a member of the nuclear receptor superfamily. It is a primary regulator of calcium and phosphate homeostasis required for skeleton and bone mineralization. Vitamin D in active form 1α,25 dihydroxyvitamin-D3 mediates its cellular functions by binding to VDR. Active VDR forms heterodimers with partner RXR (retinoid X receptor) to execute its physiological actions. HVDRR (Hereditary Vitamin D-Resistant Rickets) is a rare genetic disorder that occurs because of generalized resistance to the 1α,25(OH)2D3. HVDRR is caused by the polymorphic variations in VDR gene leading to defective intestinal calcium absorption and mineralization of newly forming bones. Using point and deletion SNPs of VDR we have studied several HVDRR-associated SNP variants for their subcellular dynamics, transcriptional functions, 'genome bookmarking', heterodimeric interactions with RXR, and receptor stability. We previously reported that VDR is a 'mitotic bookmarking factor' that remains constitutively associated with the mitotic chromatin to inherit 'transcriptional memory', however the mechanistic details remained unclear. We document that 'genome bookmarking' property by VDR is critically impaired by naturally occurring HVDRR-associated point and deletion variants found in patients. Furthermore, these HVDRR-associated SNP variants of VDR were found to be compromised in transcriptional function, nuclear translocation, protein stability and intermolecular interactions with its heterodimeric partner RXR. Intriguingly, majority of these disease-allied functional defects failed to be rescued by RXR. Our findings suggest that the HVDRR-associated SNP variations influence the normal functioning of the receptor, and this derived understanding may help in the management of disease with precisely designed small molecule modulators.

Integrative multiomics and in silico analysis revealed the role of ARHGEF1 and its screened antagonist in mild and severe COVID-19 patients.

COVID-19 is a sneaking deadly disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The rapid increase in the number of infected patients worldwide enhances the exigency for medicines. However, precise therapeutic drugs are not available for COVID-19; thus, exhaustive research is critically required to unscramble the pathogenic tools and probable therapeutic targets for the development of effective therapy. This study utilizes a chemogenomics strategy, including computational tools for the identification of viral-associated differentially expressed genes (DEGs), and molecular docking of potential chemical compounds available in antiviral, anticancer, and natural product-based libraries against these DEGs. We scrutinized the messenger RNA expression profile of SARS-CoV-2 patients, publicly available on the National Center for Biotechnology Information-Gene Expression Omnibus database, stratified them into different groups based on the severity of infection, superseded by identification of overlapping mild and severe infectious (MSI)-DEGs. The profoundly expressed MSI-DEGs were then subjected to trait-linked weighted co-expression network construction and hub module detection. The hub module MSI-DEGs were then exposed to enrichment (gene ontology + pathway) and protein-protein interaction network analyses where Rho guanine nucleotide exchange factor 1 (ARHGEF1) gene conjectured in all groups and could be a probable target of therapy. Finally, we used the molecular docking and molecular dynamics method to identify inherent hits against the ARHGEF1 gene from antiviral, anticancer, and natural product-based libraries. Although the study has an identified significant association of the ARHGEF1 gene in COVID19; and probable compounds targeting it, using in silico methods, these targets need to be validated by both in vitro and in vivo methods to effectively determine their therapeutic efficacy against the devastating virus.

Spatio-temporal dynamics of intra-host variability in SARS-CoV-2 genomes.

During the course of the COVID-19 pandemic, large-scale genome sequencing of SARS-CoV-2 has been useful in tracking its spread and in identifying variants of concern (VOC). Viral and host factors could contribute to variability within a host that can be captured in next-generation sequencing reads as intra-host single nucleotide variations (iSNVs). Analysing 1347 samples collected till June 2020, we recorded 16 410 iSNV sites throughout the SARS-CoV-2 genome. We found ∼42% of the iSNV sites to be reported as SNVs by 30 September 2020 in consensus sequences submitted to GISAID, which increased to ∼80% by 30th June 2021. Following this, analysis of another set of 1774 samples sequenced in India between November 2020 and May 2021 revealed that majority of the Delta (B.1.617.2) and Kappa (B.1.617.1) lineage-defining variations appeared as iSNVs before getting fixed in the population. Besides, mutations in RdRp as well as RNA-editing by APOBEC and ADAR deaminases seem to contribute to the differential prevalence of iSNVs in hosts. We also observe hyper-variability at functionally critical residues in Spike protein that could alter the antigenicity and may contribute to immune escape. Thus, tracking and functional annotation of iSNVs in ongoing genome surveillance programs could be important for early identification of potential variants of concern and actionable interventions.

Association of PARK-2 Non-synonyms Polymorphisms and Their In Silico Validation Among North Indian Colorectal Cancer Patients.

PURPOSE: PARK2 is a potential tumour suppressor gene and its genetic alterations (regionic loss) are common across many human cancers. The association of PARK2 germline variations (SNPs) with Parkinson's has been shown, but their association in development and progression of cancer remains elusive. The aim of this study was to identify association of PARK2 polymorphisms (rs1801474, rs1801334) with colorectal cancer in a case control study design. METHODS: This case control study included a total of 650 genetically unrelated subjects comprising 300 colorectal cancer cases and 350 healthy controls belonging to North Indian. Both SNPs were analyzed using the PCR-RFLP assay. Statistical analysis for describing risk and association was performed using SPSS-17 software. Structural deviations due to non- synonymous substitutions (S167N and D394N) were analyzed using MD simulations. RESULTS: The genotype distributions of both the SNPs were in Hardy-Weinberg equilibrium. For both the polymorphisms, the allelic model showed statistically significant risk with OR ~ 1.3. Many of the associations remained significant even after Bonferroni correction (P < 0.00125). The result suggested that both S167N and D394N were deviated from wild type and structures and were stable after 5 ns. The average value of RMSD for backbone atoms was calculated from 5 to 10 ns molecular dynamics simulation data. CONCLUSION: In conclusion, our study revealed a significant association of PARK2 SNPs with colorectal cancer as well as their relations with other clinical parameters highlighting their contribution towards colorectal cancer susceptibility in North Indian population.

Genomic characterization and epidemiology of an emerging SARS-CoV-2 variant in Delhi, India.

Delhi, the national capital of India, experienced multiple severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) outbreaks in 2020 and reached population seropositivity of >50% by 2021. During April 2021, the city became overwhelmed by COVID-19 cases and fatalities, as a new variant, B.1.617.2 (Delta), replaced B.1.1.7 (Alpha). A Bayesian model explains the growth advantage of Delta through a combination of increased transmissibility and reduced sensitivity to immune responses generated against earlier variants (median estimates: 1.5-fold greater transmissibility and 20% reduction in sensitivity). Seropositivity of an employee and family cohort increased from 42% to 87.5% between March and July 2021, with 27% reinfections, as judged by increased antibody concentration after a previous decline. The likely high transmissibility and partial evasion of immunity by the Delta variant contributed to an overwhelming surge in Delhi.

SARS-CoV-2 B.1.617.2 Delta variant replication and immune evasion.

The B.1.617.2 (Delta) variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was first identified in the state of Maharashtra in late 2020 and spread throughout India, outcompeting pre-existing lineages including B.1.617.1 (Kappa) and B.1.1.7 (Alpha)1. In vitro, B.1.617.2 is sixfold less sensitive to serum neutralizing antibodies from recovered individuals, and eightfold less sensitive to vaccine-elicited antibodies, compared with wild-type Wuhan-1 bearing D614G. Serum neutralizing titres against B.1.617.2 were lower in ChAdOx1 vaccinees than in BNT162b2 vaccinees. B.1.617.2 spike pseudotyped viruses exhibited compromised sensitivity to monoclonal antibodies to the receptor-binding domain and the amino-terminal domain. B.1.617.2 demonstrated higher replication efficiency than B.1.1.7 in both airway organoid and human airway epithelial systems, associated with B.1.617.2 spike being in a predominantly cleaved state compared with B.1.1.7 spike. The B.1.617.2 spike protein was able to mediate highly efficient syncytium formation that was less sensitive to inhibition by neutralizing antibody, compared with that of wild-type spike. We also observed that B.1.617.2 had higher replication and spike-mediated entry than B.1.617.1, potentially explaining the B.1.617.2 dominance. In an analysis of more than 130 SARS-CoV-2-infected health care workers across three centres in India during a period of mixed lineage circulation, we observed reduced ChAdOx1 vaccine effectiveness against B.1.617.2 relative to non-B.1.617.2, with the caveat of possible residual confounding. Compromised vaccine efficacy against the highly fit and immune-evasive B.1.617.2 Delta variant warrants continued infection control measures in the post-vaccination era.

Geminivirus Betasatellite-Encoded ßC1 Protein Exhibits Novel ATP Hydrolysis Activity That Influences Its DNA-Binding Activity and Viral Pathogenesis.

Plant virus satellites are maintained by their associated helper viruses, and satellites influence viral pathogenesis. Diseases caused by geminivirus-betasatellite complexes can become epidemics and therefore have become a threat to economically important crops across the world. Here, we identified a novel molecular function of the betasatellite-encoded pathogenicity determinant ßC1. The tomato leaf curl Patna betasatellite (ToLCPaB)-encoded ßC1 protein was found to exhibit novel ATPase activity in the presence of the divalent metal ion cofactor MgCl2. Moreover, ATPase activity was confirmed to be ubiquitously displayed by ßC1 proteins encoded by diverse betasatellites. Mutational and sequence analysis showed that conserved lysine/arginine residues at positions 49/50 and 91 of ßC1 proteins are essential for their ATPase activity. Biochemical studies revealed that the DNA-binding activity of the ßC1 protein was interfered with by the binding of ATP to the protein. Mutating arginine 91 of ßC1 to alanine reduced its DNA-binding activity. The results of docking studies provided evidence for an overlap of the ATP-binding and DNA-binding regions of ßC1 and for the importance of arginine 91 for both ATP-binding and DNA-binding activities. A mutant betasatellite with a specifically ßC1-ATPase dominant negative mutation was found to induce symptoms on Nicotiana benthamiana plants similar to those induced by wild-type betasatellite infection. The ATPase function of ßC1 was found to be negatively associated with geminivirus-betasatellite DNA accumulation, despite the positive influence of this ATPase function on the accumulation of replication-associated protein (Rep) and ßC1 transcripts. IMPORTANCE Most satellites influence the pathogenesis of their helper viruses. Here, we characterized the novel molecular function of ßC1, a nonstructural pathogenicity determinant protein encoded by a betasatellite. We demonstrated the display of ATPase activity by this ßC1 protein. Additionally, we confirmed the ubiquitous display of ATPase activity by ßC1 proteins encoded by diverse betasatellites. The lysine/arginine residues conserved at positions 49 and 91 of ßC1 were found to be crucial for its ATPase function. DNA-binding activity of ßC1 was found to be reduced in the presence of ATP. Inhibition of ATPase activity of ßC1 in the presence of an excess concentration of cold ATP, GTP, CTP, or UTP suggested that the purified ßC1 can also hydrolyze other cellular nucleoside triphosphates (NTPs) besides ATP in vitro. These results established the importance of the ATPase and DNA-binding activities of the ßC1 protein in regulating geminivirus-betasatellite DNA accumulation in the infected plant cell.

Phylogenetic Relationships and Potential Functional Attributes of the Genus Parapedobacter: A Member of Family Sphingobacteriaceae.

The genus Parapedobacter was established to describe a novel genus within the family Sphingobacteriaceae and derives its name from Pedobacter, with which it is shown to be evolutionarily related. Despite this, Parapedobacter and Pedobacter do not share very high 16S rRNA gene sequence similarities. Therefore, we hypothesized whether these substantial differences at the 16S rRNA gene level depict the true phylogeny or that these genomes have actually diverged. Thus, we performed genomic analysis of the four available genomes of Parapedobacter to better understand their phylogenomic position within family Sphingobacteriaceae. Our results demonstrated that Parapedobacter is more closely related to species of Olivibacter, as opposed to the genus Pedobacter. Further, we identified a significant class of enzymes called pectinases with potential industrial applications within the genomes of Parapedobacter luteus DSM 22899T and Parapedobacter composti DSM 22900T. These enzymes, specifically pectinesterases and pectate lyases, are presumed to have largely different catalytic activities based on very low sequence similarities to already known enzymes and thus may be exploited for industrial applications. We also determined the complete Bacteroides aerotolerance (Bat) operon (batA, batB, batC, batD, batE, hypothetical protein, moxR, and pa3071) within the genome of Parapedobacter indicus RK1T. This expands the definition of genus Parapedobacter to containing members that are able to tolerate oxygen stress using encoded oxidative stress responsive systems. By conducting a signal propagation network analysis, we determined that BatD, BatE, and hypothetical proteins are the major controlling hubs that drive the expression of Bat operon. As a key metabolic difference, we also annotated the complete iol operon within the P. indicus RK1T genome for utilization of all three stereoisomers of inositol, namely myo-inositol, scyllo-inositol, and 1D-chiro-inositol, which are abundant sources of organic phosphate found in soils. The results suggest that the genus Parapedobacter holds promising applications owing to its environmentally relevant genomic adaptations, which may be exploited in the future.

Hamiltonian energy as an efficient approach to identify the significant key regulators in biological networks.

The topological characteristics of biological networks enable us to identify the key nodes in terms of modularity. However, due to a large size of the biological networks with many hubs and functional modules across intertwined layers within the network, it often becomes difficult to accomplish the task of identifying potential key regulators. We use for the first time a generalized formalism of Hamiltonian Energy (HE) with a recursive approach. The concept, when applied to the Apoptosis Regulatory Gene Network (ARGN), helped us identify 11 Motif hubs (MHs), which influenced the network up to motif levels. The approach adopted allowed to classify MHs into 5 significant motif hubs (S-MHs) and 6 non-significant motif hubs (NS-MHs). The significant motif hubs had a higher HE value and were considered as high-active key regulators; while the non-significant motif hubs had a relatively lower HE value and were considered as low-active key regulators, in network control mechanism. Further, we compared the results of the HE analyses with the topological characterization, after subjecting to the three conditions independently: (i) removing all MHs, (ii) removing only S-MHs, and (iii) removing only NS-MHs from the ARGN. This procedure allowed us to cross-validate the role of 5 S-MHs, NFk-B1, BRCA1, CEBPB, AR, and POU2F1 as the potential key regulators. The changes in HE calculations further showed that the removal of 5 S-MHs could cause perturbation at all levels of the network, a feature not discernible by topological analysis alone.

Phosphatidic acid homeostasis regulated by a type-2 phosphatidic acid phosphatase represents a novel druggable target in malaria intervention.

Type-2 phosphatidic acid phosphatase (PAP2) a member of PAP2 superfamily mediates the conversion of phosphatidic acid (PA) to diacylglycerol (DAG) and thus plays a pivotal role in numerous cellular signaling processes in diverse organisms. An elevated level of intracellular PA is detrimental for the cell and induces cell death. In this study we identified and characterized a PAP2 homologue in Plasmodium falciparum, PfPAP2 and further elucidated its significance in regulation of PA homeostasis in parasite life cycle. PfPAP2 is expressed in the blood stage and harbors the canonical acid phosphatase domain (APD) with signature motifs. PfPAP2 catalyzes the dephosphorylation of PA to produce DAG and inorganic phosphate (Pi). Propranolol, a generic inhibitor of PAP2, inhibited the phosphatase activity of PfPAP2 by binding to the active site of APD domain as evident by in silico docking and confirmed by surface plasmon resonance (SPR) analysis. Inhibition of native PfPAP2 by propranolol led to rise in intracellular PA mediating disruption of intracellular PA homeostasis in parasites. The propranolol mediated inhibition of PfPAP2 directed early secretion of a micronemal Perforin like Protein, PfPLP1 leading to untimely permeabilization and host cell egress. The merozoites following premature egress were non-invasive and were attenuated to invade erythrocytes and cannot continue next cycle growth. This study demonstrates that disruption of PA homeostasis can cause growth retardation in malaria parasites, and thus its master regulator, PfPAP2, can serve as a very good molecular target for antimalarial chemotherapeutic interventions.